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The aging dynamics of colloidal suspensions of Laponite, a synthetic clay, is investigated using dy- 
namic light stattering (DLS) and viscometry after a quench into the glassy phase. DLS allows to 
follow the diffusion of Laponite particles and reveals that there are two modes of relaxation. The 
fast mode corresponds to a rapid diffusion of particles within "cages" formed by the neighboring 
particles. The slow mode corresponds to escape from the cages: its average relaxation time increases 
exponentially fast with the age of the glass. In addition, the slow mode has a broad distribution of 
relaxation times, its distribution becoming larger as the system ages. Measuring the concomitant 
increase of viscosity as the system ages, we can relate the slowing down of the particle dynamics to 
the viscosity. 

PACS numbers: 71.55.Jv, 64.70.Pf, 61.20.Lc 



I. INTRODUCTION 

Glassy systems are characterized by an equilibration 
time that is much longer than any laboratory time scale: 
glasses are out of equilibrium systems. One of the inte- 
resting consequences is that, in any case for short times 
after a quench into the glassy phase, the properties of 
a glass may still evolve on an observable time scale: we 
say that the system ages. The typical relaxation time 
of the glass does not become infinite immediately, but 
grows with the waiting time t Wl i.e., the time expired 
after the quench into the glassy phase was done. Under- 
standing the aging processes is of great interest as the 
glass transition is still ill-understood. 

Recent progress in mode-coupling theory has allowed 
for a first detailed description of the aging process [Q. 
The aging dynamics of the system is described in term 
of the evolution of the correlation and response functions 
with the waiting time. In most systems such as structural 
glasses, unfortunately correlation and response functions 
are not easy to obtain experimentally. Testing the appli- 
cability of the new mode coupling theory has therefore 
been limited to specific theoretical models and simula- 
tions (Lennard- Jones glass) studies |0,p| . 

Colloidal glasses can be looked upon as model glasses 
compared to structural glasses where there are two cou- 
pled control parameters, density and temperature. In 
colloidal glasses, the volume fraction of the particles is 
one of the control parameters that can be varied inde- 
pendently from the others. Because of the experimental 
difficulties in studying structural glasses, colloidal glasses 
have been studied extensively. In fact, it is probably fair 
to say that the strongest evidence for the applicability of 
classical mode coupling theory (that applies only above 
the glass transition temperature and therefore cannot de- 
scribe aging) comes from colloidal systems Hf|. 

In order to see whether colloidal glasses can also be 



used to investigate aging processes, we study colloidal 
suspensions of Laponite, a synthetic clay composed of 
monodisperse disc-shaped particles. These colloidal sus- 
pensions form very viscous phases at low particle con- 
centrations |^||. The first evidence of aging in colloidal 
glasses of Laponite was reported in 0. In this paper, 
we present quantitative results of aging dynamics in a 
glassy suspension of Laponite using dynamic light scat- 
tering measurements and viscometry. 



II. PREPARATION AND VISCOMETRY 

The experiments are performed with aqueous suspen- 
sions of Laponite RD, a hectorite synthetic clay provided 
by Laporte Ltd. The particles are colloidal disks of 25 
nm diameter and 1 nm thickness § , with a negative sur- 
face charge on both faces. Laponite powder was mixed 
with ultra-pure water at pH =10 obtained by addition of 
NaOH, providing chemically stable particles |9|JlO|]. The 
suspension was stirred vigorously during 15 minutes and 
then filtered through a Millipore Millex-AA 0.8 /im fil- 
ter unit. This preparation procedure allows us to obtain 
a reproducible initial liquid state. Evaporation of wa- 
ter or CO2 contamination of the sample was completely 
avoided by covering the sample with a layer of vaseline 
oil. The aging time t w — is defined as the moment the 
suspension passes through the filter. 

Suspensions prepared in this way are optically trans- 
parent and are initially liquid. Within a time varying 
from a few minutes to a few hours for the different con- 
centrations considered here (2% wt to 5% wt), a three 
order of magnitude increase in viscosity is observed : the 
suspension does not flow anymore when tumbling the re- 
cipient. Since the physical properties of the dispersion 
depend on the time after preparation t w , the conclusion 
is that the samples age. That the origin of this aging 



1 



behavior is due to the glassy dynamics was shown by 
Bonn et al. [p]Jll|. They showed that Laponite suspen- 
sions form glasses and not gels, as had been proposed 
before Jl(],[lj|]. Bonn et al. |ll| studied the structure of 
Laponite suspensions using static light scattering. Con- 
trary to previous observations, they found no evidence 
for a fractal- like organization of the particles, provided 
the suspensions were filtered. They consequently showed 
that the large viscosity increase is not due to the forma- 
tion of a fractal network. In addition, it has been shown 
recently that Laponite powder dissolves quickly into both 
individual particles and aggregates of a few particles dis- 
persing more slowly jl^] . At the final state of dispersion 
(about 5 hours), the system contains aggregates respon- 
sible for the strong scattering at small scattering vectors. 
In our experiments, aggregates of particles are broken up 
by the filtration. 

This colloidal glass is obtained for very low volume 
fractions $ ~ 0.01 compared to those for usual spherical 
colloids, for which glasses are obtained above $ ~ 0.5 [|J. 
To account for this difference, it was proposed that the 
Laponite suspensions form so-called Wigner glasses, very 
low-density glasses whose formation is due to the exis- 
tence of long-range Coulombic repulsions [p^|l6| . These 
repulsions originate from the strong surface charges at 
the faces of the colloidal disks, making that the effective 
volume fraction is very high. Recent experiments have 
shown that the location of the "glass transition" line in 
the (volume fraction / electrolyte concentration) phase 
diagram is consistent with this assumption []l7f . 

The Laponite suspensions age on time scales that de- 
pend on the particle concentration. Oscillatory shear ex- 
periments were performed on a controlled stress (Reo- 
logica StressTech) rheometer using a Couette geometry 
with a 1 mm gap. The sample was exposed to a sinu- 
soidal strain 7 = 7osinaj/j, with 70 = 0.01 and u> = 1 
s _1 , and allowing for measurement of the visco-elastic 
response. The storage modulus G' and the loss modulus 
G" were determined in the linear visco-elastic regime. 
At t w — 0, the loss and storage moduli are roughly of 
the same order of magnitude with G" > G' . As the 
system ages, the storage modulus increases more rapidly 
than the loss modulus, becoming easily more than two 
orders of magnitude larger within a typical experimental 
time scale. Figure |l| shows the evolution of the com- 
plex viscosity 77*, as a function of the aging time t w , for 
suspensions with different concentrations. The complex 
viscosity rf w as calculate d from the visco-elastic mod- 
uli, as 77* = \J G' 2 + G" 2 /uj. The complex viscosity of 
the suspension with 2.5% wt increases by three orders 
of magnitude on a time scale of 100 minutes while the 
same increase of viscosity arises on about 10 minutes for 
a suspension with 3.5% wt. The aging dynamics of the 
suspensions can therefore be varied by changing the par- 
ticle concentration. 



III. DYNAMIC LIGHT SCATTERING 

The dynamics of formation of the colloidal glass was 
investigated using dynamic light scattering (DLS). The 
light scattering set-up is as follows: a He-Ne laser beam 
(A = 632.8 nm) is focused on a cylindrical sample in 
an index-matching bath of Toluene (n — 1.50). The 
scattered intensity transmitted through the cell is de- 
tected by an optical fiber coupled to an avalanche pho- 
todiode under a scattering angle 9, defined as the angle 
between the transmitted and the scattered beam. The 
signal is then analysed by a ALV-5000 logarithmic cor- 
relator, which directly calculates the normalized inten- 
sity autocorrelation function over nine decades in time 
2 (q,*) - 1 =< J(q, 0)I(q, t) > / < J(q,0) > 2 , where t 
is the delay time. The modulus of the scattering wave 
vector, defined as q = (Ann/ \)sin(9/2), with n the re- 
fractive index of the suspension (n = 1.33) was between 
9.10 6 m _1 < q < 3.10 7 m -1 , as the scattering angle was 
varied between 40° and 150°. Its inverse, 2w/q, deter- 
mines the length scale probed in a DLS experiment. 

In our experiments, the intensity correlation function 
was recorded while the system is aging. In such an ex- 
periment, the acquisition time should allow for a good 
average of the auto-correlation function data without a 
significant aging of the system occuring during the exper- 
iment. Because of that, experiments were performed with 
suspensions of sufficiently low concentration so that the 
aging is slow compared to the time necessary to obtain 
the correlation functions. In the following, we present 
results obtained with a 2.5% wt suspension. The acquisi- 
tion time was chosen to be 60 seconds, allowing for a good 
average of the whole curve. We experimentally checked 
that aging was negligible during this acquisition time for 
the 2.5% wt suspension by recording the auto-correlation 
functions just before and just after the measurement. 

Auto-correlation functions taken for various aging 
times t w , under a scattering angle 6 = 90°, are shown in 
Figure |[ Mainly two relaxations can be observed. The 
first one, observed for short delay times t, is relatively 
fast and appears to be independent on the aging time. 
The second relaxation, observed for long delay times t, 
depends strongly on the waiting time. 

In order to describe the two processes quantitatively, 
92 (q, t) — 1, was fitted by a sum of an exponential and a 
stretched exponential function as: 

,g 2 (q, t) - 1 = Aexp -(t/n) + (1 - A) cxp -(t/r 2 ) a 

The stretched exponential is used since it has been found 
empirically that it provides a good description of the slow 
relaxation processes encountered in glassy systems. The 
fits corresponding to different aging times t w are shown 
in Figure g and describe the correlation functions very 
well, for all aging times. 

The first term of the fit function corresponds to a fast 
relaxation. The parameter t± was in fact determined in- 
dependently by a linear fit of ln((? 2 (q, t) — 1) for short 
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times t, in order to constrain the fitting procedure. This 
relaxation time was found to be independent on the aging 
time t w . Furthermore, studying the angular dependence 
of the scattering, it is found that the inverse of t\ varies 
as q 2 as shown in Figure [| This shows that, for short 
times, the Laponite particles undergo 'normal' Brown- 
ian motion: this q-dependence is the same as that found 
for very dilute colloidal dispersions [^9|; for our system, 
we find a collective diffusion coefficient D — 9.5 10~ 12 
m 2 s _1 , of the same order of magnitude as the one mea- 
sured in a very dilute solution [Q. On the other hand, 
the relaxation time T2, that characterizes the aging, was 
found to increase exponentially fast as a function of t w as 
7"2 = To ex P (tw/to), where To ~ 0.1 ms and to ~ 10 s 
(Figure ||; see also the discussion below). As the total 
decay time of the correlation function is a measure for 
the time a particle needs to "forget" its initial position, 
this shows that very rapidly, the aging freezes in a certain 
degrees of freedom of the system. 

As was the case for t\ , the relaxation time t 2 was found 
also to scale with 1/q 2 , a feature again reminiscent of 
classical (Fickian) diffusion, although the stretched ex- 
ponential shows that there are strong hydrodynamic in- 
teractions between the particles. Finally, the stretch ex- 
ponent a was found to depend only on the aging time 
and not on the scattering wavevector q. The exponent a 
decreases roughly linearly between 1 and 0, as shown in 
Figure ||. The aging time t w for which a ~ corresponds 
to the ergodic-nonergodic transition on the laboratory 
time scale Rj. 

The observation of the scaling of t\ and t 2 with 1/q 2 
allows us to rescale the autocorrelation functions for dif- 
ferent scattering angles. Rescaling the time t as q 2 a 2 t, 
where a is the radius of the Laponite particles, leads to a 
collapse of the different correlation functions (Figure ^) . 

In order to interpret these results, we determine the 
distribution of relaxation times r in the system. This we 
do using a commercially available constrained regular- 
ization method (ALV-NonLin Data Analysis p8j). The 
method consists in decomposing the intensity autocorre- 
lation function in exponential modes by directly invert- 
ing : g 2 (q,t) - 1 = ( fi™' exp (~Tt)G(T)dT) 2 , yield- 
ing the distribution function of decay times G(t), where 
G(t) = G(l/r). Figure shows the decay time distri- 
bution functions G(r) corresponding to different aging 
times t w = and t w = 400 min. One observes two 
distinct modes, the fast one corresponding to the first 
relaxation time t% , and the slow one to T2 . The intensity 
fast mode has been arbitrarily normalised to unity for 
the sake of comparison. This fast mode remains almost 
unchanged as the system ages. The slow relaxation ex- 
hibits two different features. First, the maximum of the 
distribution shifts to larger times, which corresponds to 
the increase in t 2 already observed from the direct fit of 
the correlation function. Second, the distribution of re- 
laxation times G(t) becomes wider as the system ages. 
This corresponds in fact to the decrease in the stretch ex- 



ponent a: the distribution of relaxation times r becomes 
wider and wider. These two combined effects make that, 
very rapidly, the longest relaxation time in the system 
exceeds the laboratory time scale of observation. The 
direct consequence of this is that, for short aging times 
t w , the correlation function decays to zero as shown in 
Figure |^. For longer aging times (t w > 10 hours for a 
2,5% wt suspension), the auto-correlation function does 
not decay to zero within the observation time scale: the 
system is no longer ergodic . 

Experimentally, the dependence of t 2 and rf on t w 
leads to the consequence that the viscosity of the suspen- 
sion is not simply proportional to t 2 as usually assumed 
in simulations |p8|-p0t (Figure ||). 



IV. DISCUSSION AND CONCLUSION 

We propose the following interpretation of our mea- 
surements. The diffusion of particles can be described 
gc-diffusion process. The first relaxation n char- 
acterizes the short-time Brownian diffusion of a particle 
in the suspending liquid. For short times i, the particle 
diffuses freely within a 'cage' formed by the surrounding 
particles. This diffusive motion consequently does not 
depend on the aging time. The second relaxation pro- 
cess, occurring for long times t, can be interpreted as the 
escape from the cages. The corresponding characteristic 
time t 2 increases rapidly with t w , indicating that it be- 
comes more and more difficult for a particle to 'escape'. 
This characteristic relaxation time was found also to scale 
as 1/q 2 , reminescent of diffusive motion. However, this 
diffusive motion of particles for long times is complex as 
it is characterized by a broad distribution of relaxation 
times, leading to the stretched exponential behavior of 
the correlation function. From this it follows immedi- 
ately that for short aging times, the system is ergodic: 
the particles reside in the 'cages' formed by surrounding- 
particles and escape from them after a characteristic time 
that depends on the aging time t w . As the system ages, 
the escape from the 'cages' becomes slower. The particles 
are subsequently constrained by the 'cages', resulting in 
an ergodicity breaking within the observation time scale. 

The comparison of our experimental results with 
predictions from MCT [EJ shows a qualitative agree- 
ment with respect to the slow relaxation (so-called a- 
relaxation) and faster relaxation (so-called /3-relaxation) 
features of the particle dynamics. The non- aging (3- 
relaxation part would correspond to the relaxation time 
Ti which is independent of the aging time t w . Only 
the slow a-relaxation part of the autocorrelation func- 
tion is predicted to depend on t w , which is consistent 
with the dependence of t 2 on t w . As in the MCT sce- 
nario, we found that the typical relaxation time t 2 in- 
creases with increasing t w , which leads to a non-ergodic 
system when the largest t 2 is comparable to the observa- 
tion time scale. In other experiments on aging [^3] p5[ , 
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the systems were found to relax on typical time scales 
that increase as t^,, in agreement with calculations from 
mean field models of aging p^] . When fx < 1, the pro- 
cess of aging is called sub-aging; this situation is encoun- 
tered in spin glasses E^], in concentrated colloidal sus- 
pension (24| or polymer glasses |2jJ. Recently, diffusing 
wave spectroscopy (DWS) measurements of the motion of 
the tracer particles in glassy Laponite suspensions have 
also been performed ^7j. The measurements show that 
the slow collective relaxation time scales linearly with the 
waiting time t w (so that fj, = 1) called full aging. Inter- 
estingly, in our experiments, the typical relaxation time 
T2 is found to grow exponentially fast with t w . The log- 
log inset plot of Figure || shows that indeed our data are 
inconsistent with a power law t-i oc tfo. In our experi- 
ments, the typical decay time Ti of the slow relaxation 
is of the order of 1 to 10 ms and is consequently seven 
orders of magnitude smaller than the aging time t w . In 
the experiments cited above, the aging processes were 
observed in a different regime where the aging times t w 
and the decay times r of the slow relaxation are of the 
same order of magnitude. 

In conclusion, the collective diffusion process of parti- 



cles can be interpreted ge-diffusion process. For 

short aging times, the system is ergodic. The particles 
reside in the cages formed by surrounding particles and 
escape from them after a characteristic time which is ag- 
ing time dependent. As the system ages, the escape from 
the cages become slower, resulting in a strong ergodic- 
ity breaking. In many systems, such as spin-glasses and 
concentrated colloidal supensions, the typical relaxation 
time is found to behave as t$, with \i < 1. Interestingly, 
in our experiments, the characteristic time of the slow 
relaxation process is found to grow even faster, i.e. ex- 
ponentially fast with the aging time. In addition, we 
find that this dependence on t w results in a non-trivial 
relation between the viscosity and the relaxation time. 
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FIG. 1. Complex viscosity as a function of the aging time t w for Laponite suspensions at 2.5% wt (circles) and 3.5% wt 
(squares). The visco-elastic moduli are measured on a Reologica Stress- Tech rheometer at a pulsation w — 1 s^ 1 and with an 
imposed strain 70 = 0.01. 
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FIG. 2. Autocorrelation functions of the scattered intensity by a 2.5% wt suspension and corresponding fits for aging times 
t w — 0, 200, 400, 600 min (from left to right) at a scattering angle 9 — 90°. The autocorrelation functions and their corresponding 
fits are also represented on a log-log plot in inset. 
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FIG. 3. Inverse of the characteristic time n of the fast relaxation as a function of the modulus of the wave vector q, for a 
2.5% wt Laponite suspension. 
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FIG. 4. Characteristic time T2 of the slow relaxation as a function of the aging time for 2.5% wt (squares) and 2.8% wt 
(circles) Laponite suspensions from auto-correlation functions taken under a scattering angle 9 = 90°. 7"2 (iw ) is also represented 
on a log-log plot in inset. 
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FIG. 5. Stretch exponent a as a function of the aging time for 2.5% wt (diamond) and 2.8% wt (open squares) Laponite 
suspensions. These fit parameters correspond to auto-correlation functions taken under a scattering angle 6 — 90°. 
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FIG. 6. Scaling of the autocorrelation functions taken under different scattering angles 8 = 40° (open squares), 
(open circles) and 8 = 150° (filled circles), for an aging time t w — 300 min. a is the radius of the particles M 
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FIG. 7. Relaxation time distribution functions corresponding to aging times t w = (full line) and t w = 400 minutes (dashed 
line) for a 2.5% wt suspension. The corresponding auto-correlation functions are taken under a scattering angle 6 = 90°. In 
order to be able to compare the functions, the peak intensities for the small-scale Brownian motion were normalised to unity. 
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FIG. 8. Complex viscosity as a function of T2 for a a 2, 5% wt Laponite suspension for a pulsation w = 1 s _1 . 
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